Experimental study of temperature fluctuations in forced stably stratified turbulent 

flows 
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We study experimentally temperature fluctuations in stably stratified forced turbulence in air 
flow. In the experiments with an imposed vertical temperature gradient, the turbulence is produced 
by two oscillating grids located nearby the side walls of the chamber. Particle Image Velocimetry 
is used to determine the turbulent and mean velocity fields, and a specially designed temperature 
probe with sensitive thermocouples is employed to measure the temperature field. We found that the 
ratio \{£ x V X T) 2 + (£ y V y T) 2 + (£ z V Z T) 2 ] /(9 2 ) is nearly constant, is independent of the frequency of 
the grid oscillations and has the same magnitude for both, stably and unstably stratified turbulent 
flows, where ti are the integral scales of turbulence along x, y, z directions, T and 6 are the mean and 
turbulent parts of the fluid temperature. The theoretical predictions based on the budget equations 
for turbulent kinetic energy, turbulent potential energy (oc {0 2 )) and turbulent heat flux, are in 
a good agreement with the experimental results. Detailed comparison with experimental results 
obtained previously in unstably stratified forced turbulence is performed. 

PACS numbers: 47.27.te, 47.27.-i 



I. INTRODUCTION 

In the last two decades, the theory of stably strati- 
fied flows undergoes essential revision. Since the classical 
papers by Kolmogorov, Obukhov and Hciscnberg, practi- 
cally used turbulence closures for the neutrally and stably 
stratified flows describe the energetics of turbulence us- 
ing the budget equation for the turbulent kinetic energy 
(TKE) in combination with the Kolmogorovs hypotheses 
for the dissipation rate ~ E^ 2 /£ , the eddy viscosity (or 
the eddy conductivity and eddy diffusivity) proportional 

to ~ E 1 I [ 2 £ (see Refi), where Ek is the turbulent kinetic 
energy and t is the mixing length scale. The straight- 
forward application of this approach for stably stratified 
shear flows led to the turbulence cut off at Richardson 
numbers, Ri = (N/S) 2 , exceeding some critical value, 
assumed to be close to the conventional linear instability 
threshold from 1/4 to 1 (see, e.g., Refs. 2,3 ). The lat- 
ter assertion, however, contradicts to experimental evi- 
dence and experience from numerical modelling (see, e.g., 
Refs4^— ). Here N is the Brunt- Vaisala frequency deter- 
mined by the vertical derivative of the mean temperature 
(or the mean potential temperature), i.e., iV 2 = /3V Z T, 
P = g/T* is the buoyancy parameter, g = 9.81 m/s is the 
acceleration due to gravity, T* is a reference value of the 
absolute mean temperature T and S is the mean shear. 

Over decades in meteorology, this difficulty was over- 
come heuristically importing empirical Ri-dependent co- 
efficients in the expressions for the eddy viscosity and 
eddy conductivity. Recently an insight into this long- 
standing problem (since Richardson 9 ) has been gained 
through more rigorous analysis of the turbulent ener- 
getics involving additional budget equation for the tur- 
bulent potential energy (TPE) conceptually similar to 
the Lorenzs available potential energy^, and accounting 
for the energy exchange between TKE and TPE (see, 



Refsjii— ). This analysis uses the conservation law for 
the total turbulent energy (TTE = TKE+TPE) and 
opens new prospects toward developing consistent and 
practically useful turbulent closures based on a minimal 
set of equations. This approach results in the asymp- 
totically linear Ri-dependence of the turbulent Prandtl 
number and removes the puzzling, almost the-century 
old problem of the unrealistic turbulence cut off (imply- 
ing the existence of a critical Richardson number). De- 
tailed discussions of the state of the art in the turbulence 
closure problem for stably stratified flows can be found in 
Refs* 1 ^— . These new ideas should be comprehensively 
investigated and validated using laboratory experiments 
and numerical simulations in different set-ups. 

The goal of this study is to conduct a comprehensive 
experimental investigation of heat transport in tempera- 
ture stratified forced turbulence. In the experiments tur- 
bulence is produced by the two oscillating grids located 
nearby the side walls of the chamber. We use Particle 
Image Velocimetry to determine the velocity field, and a 
specially designed temperature probe with sensitive ther- 
mocouples is employed to measure the temperature field. 
Similar experimental set-up and data processing proce- 
dure were used previously in the experimental study of 
different aspects of turbulent convection (see Refsi 17 i 18 ) 
and in Refs»i2r— for investigating a phenomenon of tur- 
bulent thermal diffusion! 24 : 25 Comprehensive investiga- 
tion of turbulent structures, mean temperature distribu- 
tions, velocity and temperature fluctuations can eluci- 
date a complicated physics related to particle clustering 
and formation of large-scale inhomogeneities in particle 
spatial distributions in stably stratified turbulent flows. 

In the present study we perform a detailed compar- 
ison with experimental results obtained recently in un- 
stably stratified forced turbulence 1 -^, whereby transition 
phenomena caused by the external forcing (i.e., transition 
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from Rayleigh-Benard convection with the large-scale cir- 
culation (LSC) to the limiting regime of unstably strati- 
fied turbulent flow without LSC where the temperature 
field behaves like a passive scalar) have been studied. In 
particular, when the frequency of the grid oscillations 
is larger than a certain value, the large-scale circula- 
tion in turbulent convection is destroyed, and the de- 
struction of the LSC is accompanied by a strong change 
of the mean temperature distribution. However, in all 
regimes of the unstably stratified turbulent flow the ratio 
[(4V,T) 2 + [lyWyTf + (£ z V z T) 2 ]/(9 2 ) varies slightly 
(even in the range of parameters whereby the behaviour 
of the temperature field is different from that of the pas- 
sive scalar) jA£ 

This paper is organized as follows. The theoretical pre- 
dictions are given in Section II. Section III describes the 
experimental set-up and instrumentation. The results of 
laboratory study of the stably stratified turbulent flow 
and comparison with the theoretical predictions are de- 
scribed in Section IV. Finally, conclusions are drawn in 
Section V. 



II. THEORETICAL PREDICTIONS 

In our theoretical analysis we use three budget equa- 
tion for the turbulent kinetic energy Ek = (u 2 )/2, for 
the temperature fluctuations Eg = {0 2 )/2 and for the 
turbulent heat flux Fi = {ui9): 



tion for the turbulent potential energy E p = [fi 2 /N 2 ) Eg: 
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(see, e.g., RefsJ^ 2 ^), w h ere D/Dt = d/dt + U-V, u 
are the fluctuations of the fluid velocity, 6 are the temper- 
ature fluctuations, U is the mean velocity, T is the mean 
temperature, p are the pressure fluctuations, ft = /3 ej, e 
is the vertical unit vector and p is the fluid density. The 
terms <frg and include the third-order moments. 
In particular, <I>fc = p _1 (up) + (l/2)(uu 2 ) determines the 
flux of Ek, 4>0 = {u9 2 )/2 determines the flux of Eg and 
— (uiUjO) + Sijp^ 1 (9p)/2 determines the flux of 
F. The term (u-f/) in Eq. ([T]) determines the production 
of turbulence caused by the grid oscillations and s k is 
the dissipation rate of the turbulent kinetic energy. The 
term eg in Eq. ^ determines the dissipation rate of Eg , 

while the term in Eq. ([3]) is the dissipation rate of 
the turbulent heat flux. 

By means of Eq. @ we arrive at the evolutionary equa- 
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(see, e.g., Refs^i*), where N 2 = /3V Z T, * p = 
(/3 2 /7V 2 )* , P p = -(p 2 /N 2 )(F h -V)T is the source (or 
sink) of the turbulent potential energy caused by the 
horizontal turbulent heat flux F/ t = (u/j#), u/ t is the 
horizontal component of the velocity fluctuations and 
e p = (0 1 /N 2 ) eg. The buoyancy term, j3 F z , appears in 
Eqs. (fT]) and (TJJ with opposite signs and describes the en- 
ergy exchange between the turbulent kinetic energy and 
the turbulent potential energy. These two terms cancel 
in the budget equation for the total turbulent energy, 
E = Ek + E p : 

DE 

— + div * = P p - ( Ui Uj) VjUi + <u-f» - e, (5) 



where 3? = <3?/, + & p and e = Sk + £> The concept of the 
total turbulent energy is very useful in analysis of strat- 
ified turbulent flows. In particular, it allows to elucidate 
the physical mechanism for the existence of the shear pro- 
duced turbulence for arbitrary values of the Richardson 
number, and abolish the paradigm of the critical Richard- 
son number in the stably stratified atmospheric turbu- 
lence (see RefSfi2J£). 

Now we use the budget equation © for the turbu- 
lent heat flux Fi — (u; L 9). According to the estimate 
made in Refi 2 , ft (9 2 ) - p- 1 (9\7 iP ) w C e ft (9 2 ), where 
Cg is an empirical constant. In a steady-state case 
Eq. ([3]) yields the components of the turbulent heat flux 
F x = -Dl V X T, F y = -D T y \7 y T and F z = -Dj V 2 T + 
Cg C F t /3 (9 2 ), where Df — C' F t (u 2 ) with i = x,y,z 
are the turbulent temperature diffusion coefficients in 
x,y, and z directions and Cf is an empirical constant. 
Here we have taken into account that the dissipation rate 
of the turbulent heat flux is e\ F ^ — Fi/CpTo, the diago- 
nal components of the tensor (uiUj) are much larger than 
the off-diagonal components, and t x « r y ~ t z — tq. 

On the other hand, in a steady-state case Eq. @ yields 



-2r (F-V)T, 



(G) 



where we have taken into account that the dissipa- 
tion rate of Eg = (9 2 )/2 is eg = (9 2 )/2t . To esti- 
mate the dissipation rate eg, we apply the Kolmogorov- 
Obukhov hypothesis: eg « (9 2 )/2t (see, e.g., Refs^^), 
where To = (/uq is the characteristic turbulent time 
and Mo is the characteristic turbulent velocity at the 
integral turbulent scale I. Indeed, eg = -D((V#) 2 } = 
D(9 2 ) f^k 2 E g (k)dk « (6» 2 )/2t , where D is the co- 
efficient of molecular temperature diffusion, Eg{k) = 
(2/3) k^ 1 (k/ko)~ 5 / 3 is the spectrum function of the tem- 
perature fluctuations, k = t~ x , kb = i^ 1 , l\, = £/Pe 3 ^ 4 
and Pc = uqI/D is the Peclct number. The latter esti- 
mate implies that the main contribution to the dissipa- 
tion rate eg arises from very small molecular temperature 
diffusion scales th- 



Substituting the components of the turbulent heat flux 
into Eq. ([5]), we obtain the following equation: 
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const, 
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where 



[4 V*T] 2 = [{l x V x Tf + {l y V y Tf + (£ Z V Z T) 2 ] 



x[l + 2C e C F (3T 2 (\7 z T)] 



-l 



(8) 



In deriving Eq. (J7} we have neglected the terms ~ 
0[£ 3 /Lj,;£ 3 /(L^,Lu)}, where Lt and Lu are the char- 
acteristic spatial scales of the mean temperature and ve- 
locity field variations. 

Next, we derive equation for the vertical turbulent ve- 
locity in the stably stratified turbulent flow. We use the 
budget equation for the vertical turbulent kinetic energy 



E z = 



<)/2: 

DE Z 
Dt 



div 



(u,f z ) + F g + > 



(see, e.g., Refs 
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where the 
p- 1 (u z p)e+ (l/2)(uu 2 ) 



-e», (9) 

third-order mo- 
dctcrmincs the 



ment <& z = 

flux of E z , {u z f z ) is the production term, and Q z is 
the inter-component energy exchange term that accord- 
ing to the " return-to- isotropy" hypothesis™ is given by 
Q z = —C r (E z — Ek/3)/2t. Here we have neglected a 
small production term due to the weak non-uniform mean 
flow. In the steady-state Eq. © yields 



E z 



T[(u z f z )+/3F Z } 
1 + C r (l- 1/34*)' 



(10) 



where C r is an empirical constant, A z = E z /Ek is the 
vertical anisotropy parameter. Using a simple estimate 
for the vertical heat flux, F z ~ — C \J {u 2 z ) \J {6 2 ) (where 
C is the correlation coefficient), we arrive at the following 
equation for the r.m.s. of the vertical turbulent velocity 
in the stably stratified turbulent flow: 



vM) ~ ((ui) 2 )-cj. z (3y/W) 



1/2 



(11) 



where C u is an empirical constant to be determined in 
the experiment, and we have taken into account that 
the characteristic velocity for the isothermal turbulence 
{(u;) 2 )^2r{u z f z )/[l + C r (l~l/3A z )}. 



III. EXPERIMENTAL SET-UP AND 
INSTRUMENTATION 

In this section we describe the experimental set-up. 
The experiments in stably stratified turbulence have 
been conducted in rectangular chamber with dimensions 
26 x 58 x 26 cm 3 in air flow with the Prandtl number 
Pr = 0.71. The side walls of the chambers are made 
of transparent Perspex with the thickness of 1 cm. In 
the experiments turbulence is produced by two oscillat- 
ing grids. Pairs of vertically oriented grids with bars 




FIG. 1. Experimental set-up: (1) -laser light sheet; (2) - heat 
exchangers; (3) - grid driver; (4) - digital CCD camera. 



arranged in a square array (with a mesh size 5 cm) are 
attached to the right and left horizontal rods (see Fig. [IJ . 
The grids are positioned at a distance of two grid meshes 
from the chamber walls and are parallel to the side walls. 
Both grids are operated at the same amplitude of 6.1 cm, 
at a random phase and at the same frequency which is 
varied in the range from 1 Hz to 10.5 Hz. To increase the 
size of the domain with a homogeneous turbulence and 
to decrease the mean velocity, vertical partitions were 
attached to the side walls of the chamber. We use the 
following system of coordinates: z is the vertical axis, 
the y-axis is perpendicular to the grids and the xz-planc 
is parallel to the grids. The aspect ratio of the cham- 
ber H y /H z — 1.52, where H y is the size of the chamber 
along y-axis between partitions and H z is the hight of 
the chamber, respectively. 

A vertical mean temperature gradient in the turbulent 
air flow was formed by attaching two aluminium heat ex- 
changers to the bottom and top walls of the test section 
(a cooled bottom and a heated top wall of the cham- 
ber). To improve heat transfer in the boundary layers at 
the bottom and top walls we used heat exchangers with 
rectangular fins 0.3 x 0.3 x 1.5 cm 3 (see Fig. [T|) which 
allowed us to form a mean temperature gradient up to 
1.15 K/cm at a mean temperature of about 308 K when 
the frequency of the grid oscillations / > 10 Hz. The 
thickness of the aluminium plates with the fins is 2.5 cm. 
The bottom plate is a top wall of the tank with cool- 
ing water. Temperature of water circulating through the 
tank and the chiller is kept constant within 0.1 K. Cold 
water is pumped into the cooling system through two in- 
lets and flows out through two outlets located at the side 
walls of the cooling system. The top plate is attached to 
the electrical heater that provides constant and uniform 
heating. The voltage of a stable power supply applied to 
the heater varies up to 200 V. The power of the heater 
varies up to 300 W. 
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The temperature field was measured with a tem- 
perature probe equipped with twelve E-thermocouples 
(with the diameter of 0.0 f 3 cm and the sensitivity of 
~ 65 /iV/K) attached to a vertical rod with a diameter 
0.4 cm. The spacing between thermal couples along the 
rod was 2.2 cm. Each thermocouple was inserted into a 
O.f cm diameter and 4.5 cm long case. A tip of a ther- 
mocouple protruded at the length of 1.5 cm out of the 
case. The temperature in the central part of the chamber 
was measured for 2 rod positions with 2.2 cm intervals 
in the horizontal direction, i.e., at 24 locations in a flow. 
The exact position of each thermocouple was measured 
using images captured with the optical system employed 
in PIV measurements. A sequence of 600 temperature 
readings for every thermocouple at every rod position 
was recorded and processed using the developed software 
based on Lab View 7.0. 

The velocity fields were measured using a Stereoscopic 
Particle Image Velocimetry (PIV), see Refs.— ~— . In the 
experiments we used LaVision Flow Master III system. 
A double-pulsed light sheet was provided by a Nd-YAG 
laser (Continuum Surelite 2 x 170 mJ). The light sheet 
optics includes spherical and cylindrical Galilei telescopes 
with tuneable divergence and adjustable focus length. 
We used two progressive-scan 12 bit digital CCD cam- 
eras (with pixel size 6.7 /im x6.7/xm and 1280 x 1024 
pixels) with a dual-frame-technique for cross-correlation 
processing of captured images. A programmable Timing 
Unit (PC interface card) generated sequences of pulses 
to control the laser, camera and data acquisition rate. 
The software package LaVision DaVis 7 was applied to 
control all hardware components and for 32 bit image ac- 
quisition and visualization. This software package com- 
prises PIV software for calculating the flow fields using 
cross-correlation analysis. 

To obtain velocity maps in the central region of the 
flow in the cross-section parallel to the grids and perpen- 
dicular to a front view plane, we used one camera with a 
single-axis Scheimpflug adapter. The angle between the 
optical axis of the camera and the front view plane as 
well as the angle between the optical axis and the probed 
cross-section was approximately 45 degrees. The per- 
spective distortion was compensated using Stereoscopic 
PIV system calibration kit whereby the correction was 
calculated for a recorded image of a calibration plate. 
The corrections were introduced in the recorded in the 
probed cross-section images before their processing us- 
ing a cross-correlation technique for determining velocity 
fields. 

An incense smoke with sub-micron particles (p p / p ~ 
10 3 ), was used as a tracer for the PIV measurements. 
Smoke was produced by high temperature sublimation 
of solid incense grains. Analysis of smoke particles us- 
ing a microscope (Nikon, Epiphot with an amplification 
of 560) and a PM-300 portable laser particulate analyzer 
showed that these particles have an approximately spher- 
ical shape and that their mean diameter is of the order 
of 0.7/xm. The probability density function of the par- 



ticle size measured with the PM300 particulate analyzer 
was independent of the location in the flow for incense 
particle size of 0.5 — 1 /im. The maximum tracer particle 
displacement in the experiment was of the order of 1/4 
of the interrogation window. The average displacement 
of tracer particles was of the order of 2.5 pixels. The av- 
erage accuracy of the velocity measurements was of the 
order of 4% for the accuracy of the correlation peak de- 
tection in the interrogation window of the order of 0.1 
pixel (see, e.g., Refsi^Ir— ). 

We determined the mean and the r.m.s. velocities, two- 
point correlation functions and an integral scale of tur- 
bulence from the measured velocity fields. Series of 520 
pairs of images acquired with a frequency of 1 Hz, were 
stored for calculating velocity maps and for ensemble and 
spatial averaging of turbulence characteristics. The cen- 
ter of the measurement region in yz and xz-planes co- 
incides with the center of the chamber. We measured 
velocity in a flow domain 25.6 x 25.6 cm 2 with a spatial 
resolution of 1024 x 1024 pixels. This corresponds to a 
spatial resolution 250 /im / pixel. The velocity field in 
the probed region was analyzed with interrogation win- 
dows of 32 x 32 or 16 x 16 pixels. In every interrogation 
window a velocity vector was determined from which ve- 
locity maps comprising 32 x 32 or 64 x 64 vectors were 
constructed. The mean and r.m.s. velocities for every 
point of a velocity map were calculated by averaging over 
520 independent maps, and then they were averaged over 
the central flow region. 

The two-point correlation functions of the velocity field 
were determined for every point of the central part of the 
velocity map (with 16 x 16 vectors) by averaging over 
520 independent velocity maps, which yields 16 correla- 
tion functions in y and z directions. Then the two-point 
correlation function was obtained by averaging over the 
ensemble of these correlation functions. An integral scale 
of turbulence, £, was determined from the two-point cor- 
relation functions of the velocity field. In the experiments 
we evaluated the variability between the first and the last 
20 velocity maps of the series of the measured velocity 
field. Since very small variability was found, these tests 
showed that 520 velocity maps contain enough data to 
obtain reliable statistical estimates. We found also that 
the measured mean velocity field is stationary. 

The characteristic turbulence time in the experiments 
t z = 0.2 — 1 seconds that is much smaller than the time 
during which the velocity fields are measured (~ 600 s). 
The size of the probed region did not affect our results. 
The temperature difference between the top and bottom 
plates, AT, in all experiments was 50 K. 



IV. EXPERIMENTAL RESULTS AND 
COMPARISON WITH THE THEORETICAL 
PREDICTIONS 

We will start with the discussion of the experimental 
results on the temperature measurements in stably strat- 
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FIG. 2. Vertical profile of the mean temperature T(z) for 
different frequencies / of the grid oscillations for the stably 
stratified turbulent flow: / = 1 Hz (triangles), f — 2 Hz 
(squares), and / = 10.5 Hz (stars) for AT = 50 K. The mean 
temperature is measured in K and the distance in cm. 

ified turbulent flows. To avoid side effects of the grids we 
present the experimental results recorded in the central 
region of the chamber with the size of 10 x 10 x 10 cm 3 . 
The temperature was measured at 24 locations in a flow 
(the spacing between thermal couples was 2.2 cm). To 
determine the mean temperature field, T, we perform 
sliding averaging of the instantaneous temperature field 
over time that is by one order of magnitude larger than 
the characteristic turbulence time. The turbulence time 
in different experiments varies from 0.2 to 1 s depending 
on the frequency of the grid oscillations. The separation 
distance of 2.2 cm between thermal couples is sufficient 
to measure the gradients of the mean temperature. In- 
deed, the integral scale of turbulence is about 2 cm and 
the characteristic length scale of the mean temperature 
field is much larger than the integral scale of turbulence. 

Let us first analyze the frequency dependence of ver- 
tical profiles of the mean temperature T(z) (Fig. [2]). In- 
spection of Fig.[5]shows that the increase of the frequency 
/ of the grid oscillations weakly modifies the vertical pro- 
files of the mean temperature. However, the gradients of 
the mean temperature in the vertical, V Z T, and horizon- 
tal, V y T, directions are affected by the increase of the 
frequency / of the grid oscillations (Fig. [3]). 

Similarly to the unstably stratified turbulence^ we 
have observed the long-term nonlinear oscillations of 
the mean temperature with the periods which are much 
larger than the turbulent correlation time (see Fig. 0}. 
In particular, in Fig. U we show time dependencies of 
the instantaneous temperature T tot = T + 6, the long- 
term variations of mean temperature ST = T — T and 
the long-term variations of the vertical mean tempera- 
ture gradient <5(V 2 T) = V Z T - V Z T due to the nonlin- 
ear oscillations of the mean temperature. The bar here 
denotes additional time averaging. We also determined 
the long-term variations of the mean temperature gra- 
dients (5(VjT) = V,-T — ViT in other directions, where 



FIG. 3. Vertical gradients of the mean temperature V 2 T 
(squares) and r.m.s. of 5(V Z T) (diamonds) characterizing 
amplitude of the long-term nonlinear oscillations versus the 
frequency / of the grid oscillations for the stably stratified 
turbulent flow for AT = 50 K. Similar dependencies are also 
shown for horizontal gradients of the mean temperature V M T 
(snowflakes) and r.m.s. of <5(V y T) (stars). The mean temper- 
ature gradient is measured in K cm -1 and the frequency / is 
measured in Hz. 
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FIG. 4. Time dependencies of the instantaneous temperature 
T tot = T + e,the variations of mean temperature ST = T - T 
and the variations of the vertical mean temperature gradient 
<5(V Z T) = V 2 T — V Z T due to the long-term nonlinear oscil- 
lations of the mean temperature (with the period ~ 12 s). 
The bar denotes additional time averaging. These time de- 
pendencies are measured in the center of the chamber at the 
frequency / = 10.5 Hz of the grid oscillations for the stably 
stratified turbulent flow for AT = 50 K. These temperature 
characteristics are measured in K and time is measured in 
seconds. 

i — x, y, z. 

In Fig. [5] we show the results of a Fourier analysis 
of the signal ST = T — T. Inspection of Fig. [5] shows 
that there are two main maxima in the spectrum with 
the periods 12 s and 20.8 s. Other smaller maxima in 
the spectrum are at the frequencies which are multiples 
of these main frequencies or their sums and differences. 
These are typical features of nonlinear oscillations. The 
theory that explains the mechanism of these nonlinear 



6 



£«(/) 
15 

10 

5 



O.OSflz 
I 




0.10 Hz °- 13 ' 0l5 > O^tHz 




0.05 0.1 0.15 0.2 



FIG. 5. The normalized spectrum function E n (f) = 
\{8T)f\ 2 / J \(5T) f -\ 2 df of the signal ST = T(t) - T, where 

in the Fourier space (ST) / = J ST exp[— if t] dt and / is the 
frequency of the nonlinear long-term oscillations of the mean 
temperature. 
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FIG. 7. The r.m.s. turbulent velocity u rma versus frequency / 
of the grid oscillations for the stably stratified turbulent flow 
with AT = 50 K (squares) and for isothermal turbulence 
(triangles). 
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FIG. 6. The r.m.s. of temperature fluctuations 8 rms = yj (8 2 ) 
versus the frequency / of the grid oscillations for the stably 
(squares) and unstably (stars) stratified turbulent flows. 



FIG. 8. Vertical anisotropy A z — (v? z )/uxma versus frequency 
/ of the grid oscillations for the stably stratified turbulent 
flow with AT = 50 K (squares) and for isothermal turbulence 
(triangles) . 



oscillations of the mean temperature field, has not been 
developed yet. A possible mechanism for such nonlinear 
oscillations could be related to the large-scale Tollmien- 
Schlichting waves in sheared turbulent flows (see Refill) . 

In our experiments, we determined the dependence of 
the r.m.s. of the temperature fluctuations yj (0 2 ) versus 
the frequency / of the grid oscillations (see Fig. [6]), where 
6 are fluctuations of fluid temperature. The temperature 
fluctuations monotonically increase with the increase the 
frequency / of the grid oscillations (except for the higher 
frequency) due to the monotonic increase of the r.m.s. ve- 
locity fluctuations. In the case o f the unstably stratified 
turbulent flow the dependence yj (0 2 ) versus frequency is 
more involved. 

We also determined the frequency dependencies of the 
following measured turbulence parameters: the r.m.s. ve- 

1/2 

locity fluctuations, M rms = + u 2 + u 2 )] (Fig. [7}, 



the vertical anisotropy A z = (u 2 )/u rr[LS (Fig. [5]), and the 
integral scales of turbulence along horizontal y (Fig. |9]) 
and vertical z (Fig.[TU]) directions (£ y and £ z ). Except for 
the small frequencies of the grid oscillations, the horizon- 
tal integral scale of turbulence behaves similarly for both, 
the stably and unstably stratified turbulent flows, while 
the vertical integral scale of turbulence is systematically 
higher for the unstably stratified turbulent flow. On the 
other hand, the vertical anisotropy A z only slightly in- 
creases with the frequency of the grid oscillations. 

The mean, r.m.s. velocities and the two-point correla- 
tion functions of the velocity field have been calculated 
by averaging over 520 independent velocity maps, and 
then they have been averaged over the central flow re- 
gion. The integral scales of turbulence, l y and l z , have 
been determined from the normalized two-point longi- 
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FIG. 9. Horizontal integral scale of turbulence £ a versus the 
frequency / of the grid oscillations for the stably (squares) 
and unstably (stars) stratified turbulent flows with AT = 50 
K, and for isothermal turbulence (triangles). The turbulent 
length scales are measured in mm and the frequency / is 
measured in Hz. 
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FIG. 11. Reynolds number Re = lu Tma /v versus frequency 
/ of the grid oscillations for the stably stratified turbulent 
flow for AT = 50 K (squares) and for isothermal turbulence 
(triangles). 
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FIG. 10. Vertical integral scale of turbulence £ z versus the 
frequency / of the grid oscillations for the stably (squares) 
and unstably (stars) stratified turbulent flows with AT = 50 
K, and for isothermal turbulence (triangles). The turbulent 
length scales are measured in mm and the frequency / is 
measured in Hz. 



tudinal correlation functions of the velocity field, e.g., 
F v(v) = ( u v( r o)u v {r + ye v ))/(ul(r )) [and similarly 
for F z {z) after replacing in the above formula y by z], 
using the following expression: £ y = (J Q F y (y) dy)s [and 
similarly for £ z after replacing in the above formula y 
by z], where L — 10 cm is the linear size of the probed 
flow region, e y is the unit vector in the y direction, and 
(■■■)s is the additional averaging over the yz cross-section 
of the probed region. In Fig. [TT] we show the Reynolds 
number Re — lu rms jv versus frequency of the grid os- 
cillations for the stably stratified turbulent flow, where 
I = {£ 2 + £y + (.I) 1 ! 2 - The Reynolds number increases 
with the increase of the frequency of the grid oscillations 
due to increase of the production rate of turbulence. 

In Fig. [T^] we plot the non-dimensional ratio 
£* V*T/ y/ (8 2 ) [see Eq. ©] versus the frequency / of the 



FIG. 12. The non-dimensional ratio £» V,T/6 rms versus the 
frequency / of the grid oscillations for the stably (squares) 
and unstably (stars) stratified turbulent flows for AT = 50 
K, where 9 rm s = \l (0 2 ). 



grid oscillations for the stably stratified turbulent flow 
(squares) obtained in our experiments. For comparison 
in the same figure we show also this non-dimensional ra- 
tio obtained previously 17 for the unstably stratified tur- 
bulent flow (stars) . Inspection of Fig. [TJ] shows that this 
non-dimensional ratio is nearly independent of the fre- 
quency of the grid oscillations and has the same mag- 
nitude for both, stably and unstably stratified turbulent 
flows, in agreement with the theoretical predictions. Here 
we assumed that l x = £ y , Cf = 0.2 and Co = 0.83 for 
the unstably stratified turbulence, while Cg = 1.4 for the 
stably stratified turbulence. Small deviations of the ex- 
perimental results from the theoretical predictions [see 
Eq. ([7}] may be caused by a non-zero term, div <f>g and 
deviations from the steady-state due to the long-term os- 
cillations of the mean temperature. 

We also determined the frequency dependence of the 
ratio u z /u z ms (see Fig. [T5|) . where u r z ms is the r.m.s. of 
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FIG. 13. Ratios u*,/u r z ma (squares) and u%/u z (circles) ver- 
sus the frequency / of the grid oscillations. Here u™ s is the 
r.m.s. of the vertical component of velocity fluctuations in 
the unstably stratified turbulent flow, u* is the r.m.s. of the 
vertical component of velocity fluctuations in the isothermal 
turbulence, u z is the vertical component of the effective tur- 
bulent velocity, u z = [{u z ) 2 - G u l z /3 ^(6 2 }] 1/2 , that takes 
into account the attenuation of the turbulence by buoyancy. 
Here C u is an empirical constant that is about 1.81 for sta- 
bly stratified turbulence and is about 4 for unstably stratified 
turbulence. The velocity is measured in cm s _1 and the fre- 
quency / is measured in Hz. 



the vertical component of velocity fluctuations in the sta- 
bly stratified turbulent flow and u* z is the r.m.s. of the 
vertical component of velocity fluctuations in the isother- 
mal turbulence. We also determine the ratio u* z /u Zl 
where u z is the vertical component of the effective tur- 
bulent velocity, u z = [(u*) 2 - C u £ z /3 ^/W)Y /2 - Th is 
effective velocity [see Eq. (fTTj) ] that takes into account 
the destruction of the turbulence by buoyancy, is derived 
from the budget equation ^ for the vertical turbulent 
kinetic energy in Sect. II. Inspection of Fig. Q2] shows 
that the values of these ratios, u* z /u z ms and u* z /u z , are 
very close. The latter implies that the measured turbu- 
lent velocity in the stably stratified turbulent flow, u r z ms , 
is of the order of u z , in agreement with the theoretical 
predictions. 

To characterize the stably stratified flows, in Fig. Q3] 
we show the turbulent Richardson number R%t — N 2 t 2 
versus frequency / of the grid oscillations for the sta- 
bly stratified turbulent flow, where to = i/u rms is the 
characteristic turbulent time. The turbulent Richardson 
number Rix strongly decreases with the increase of the 
frequency of the grid oscillations due to the strong de- 
crease of the turbulent correlation time To with increase 
of the frequency /. For large frequencies of the grid os- 
cillations whereby Rir <C 1, the temperature field can be 
considered as a passive-scalar^ On the other hand, for 
smaller frequencies of the grid oscillations, R%t > 1, and 
the temperature field behaves as an active field. 



FIG. 14. Turbulent Richardson number Rir = N 2 t§ versus 
frequency / of the grid oscillations for the stably stratified 
turbulent flow. 



V. DISCUSSION AND CONCLUSIONS 

Temperature fluctuations in stably stratified forced 
turbulence in air flow are investigated in laboratory ex- 
periments. The stratification is caused by an imposed 
vertical temperature gradient, and the turbulence is sus- 
tained by vertical oscillating grids. It was demonstrated 
that the ratio [(£ X V X T) 2 + (^V y T) 2 + (i z \I Z T) 2 \ /(9 2 ) 
is nearly constant and is independent of the frequency of 
the grid oscillations in both, stably and unstably strati- 
fied turbulent flows. 

We found that for large frequencies of the grid oscil- 
lations the turbulent Richardson number, Rix, is small 
and the temperature field can be considered as a passive- 
scalar. Note that the passive-like scalar behaviour of 
the temperature field can be understood in the kine- 
matic sense. In particular, when the temperature fluctu- 
ations (Q 2 ) do not affect the turbulent kinetic energy, the 
temperature field can be considered as a passive scalar. 
This implies that the evolution of the temperature field 
in a given turbulent velocity field is a kinematic prob- 
lem, whereby there is no dynamic coupling between the 
temperature fluctuations, (0 2 ), and the turbulent kinetic 
energy. When the effect of the temperature fluctuations 
on the turbulent kinetic energy cannot be neglected, the 
temperature is considered as an active field. This defini- 
tion of the passive or active behavior of the temperature 
field is different from that based on the scaling behaviour 
of the temperature structure function^ 

The temperature fluctuations have been investigated 
here also theoretically using the budget equations for tur- 
bulent kinetic energy, turbulent potential energy (deter- 
mined by the temperature fluctuations) and turbulent 
heat flux. The developed theory is in a good agreement 
with the experimental results. The long-term nonlin- 
ear oscillations of the mean temperature in stably strat- 
ified turbulence have been observed for all frequencies of 
the grid oscillations similarly to the case of the unsta- 
bly stratified flow. One of the explanations of this effect 
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could be related to the large-scale Tollmien-Schlichting 
waves in sheared turbulent flows 34 , which can result in 
the nonlinear oscillations of the mean temperature field. 
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